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Abstract 

A new scaling variable is introduced in terms of which nuclear shadowing in deep- 
inelastic scattering is universal, i.e. independent of A, Q 2 and x. This variable can be 
interpreted as a measure of the number of gluons probed by the hadronic fluctuations 
of a virtual photon during their lifetime. According to recent data from the HI and 
ZEUS experiments, the gluon density in a proton rises steeply as x tends to zero. So 
nuclei which in the infinite momentum frame have a large surface density of gluons 
model at larger x values what we expect for the proton at much smaller x. Using 
experimental information on nuclear shadowing we predict the unitarity corrections 
to the proton structure function at small x and extract the bare Pomeron intercept. 
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Nuclear shadowing at small Bjorken x <C 1 is experimentally well established in deep- 
inelastic lepton scattering (DIS) [0, 0, and in Drell-Yan lepton pair production [|] . 
Theoretically shadowing in the nuclear structure functions has been predicted on the basis 
of the old fashioned parton model H |6| and more recently in QCD calculations 0, §■ In 
both approaches shadowing is the result of parton recombination. 

Shadowing, or unitarity corrections, is supposed to exist in the proton structure function 
as well. In this case, however, there is no direct way to observe experimentally shadowing 
in the proton because there is nothing to compare to. The unitarity corrections at small x 
are expected to be significant because the virtual photoabsorption cross section rises steeply 
P"T| , |TU| , steeper than hadronic cross sections at high energy. This fast growth is usually 
interpreted as being due to the Pomeron contribution to the structure function corrected 
for shadowing. No reliable way to calculate unitarity corrections has been proposed. 

Since at small x the parton clouds of different nucleons overlap ||, one expects larger 
parton densities, and consequently stronger shadowing effects, in nuclei than in protons. 
The question arises: can one use the data on nuclear shadowing at small x in order to 
improve the reliability of the unitarity corrections to the proton structure function? The 
aim of this letter is to answer this question. 

There is a widespread belief that the qluon and sea quark distribution functions expe- 
rience nearly the same nuclear shadowing, because the sea, which dominates at small x, is 
supposed to be generated through gluons. This would be true, if our electromagnetic probes 
of the gluon distribution were really hard. We argue below that this is not obvious even at 
high Q 2 . 

A photon, through its quark- ant iquark fluctuations, interacts with gluons. This process 
can also be interpreted as an electromagnetic interaction with gluonic qq fluctuations. At 
high Q 2 perturbative QCD can be applied, and it is convenient to decompose the photon 
wave function in a quark-gluon basis. Assuming that the transverse separation p of the qq 
fluctuations in a highly virtual photon is small (see, however, the discussion below) and 
that the interaction cross section vanishes as er(p) oc p 2 [I2|, one can present the total 



photoabsorption cross section in the following form |L4L 15 



Q 2 ) « ^a s (p) (p 2 ) g N (x, Q 2 ) (1) 

Here ^(x, Q 2 ) = xGn(x,Q 2 ) is the gluon distribution function. The averaging in (|l|) 
is weighted with the square of the wave function of the qq fluctuation of the photon (]1S[| , 
(p 2 ) = Jq 1 da J d 2 p\fy 2 , (p, a) |^ p 2 , where a is the fraction of the photon light-cone momentum 
carried by the quark; |\I/ 2 , (p, a)\\ oc [1 — 2ct(l — a)] e 2 Kf(ep) and e 2 = a(l — a)Q 2 + m 2 is 
an important parameter which determines the transverse separation of the qq fluctuation, 
p 2 oc 1/e 2 . It is of order 1/Q 2 except at the edges of the kinematical region, where a 
or 1 — a ~ m 2 /Q 2 , and the gg fluctuation acquires a large transverse size, p 2 ~ l/"^ 2 - 



18, 16, 17 



Of course such a component of the photon cannot be used as a probe of the gluon 
distribution. For instance, if one needs to estimate how many nucleons are in a nucleus, 
and one uses, say, a pion-nucleus inelastic interaction as a probe, one obviously will get a 
wrong result, i.e., that the number of nucleons is ~ A 2 / 3 . The source of the trouble is the 
softness of the probe, a^N) is too large. The same problem may arise in DIS. It is not 
clear what part of the observed nuclear shadowing comes from the reduction of the gluon 
density caused by recombination, and what is due to the softness of the qq probe. The latter 
corresponds to the multiple scattering shadowing mechanism. 

Although multiple scattering corrections violate the factorization property of eq. ([]]), it 
is natural to assume that both shadowing mechanisms solely depend upon the number of 
gluons involved in the interaction, rather than upon x, Q 2 or the nuclear mass number A. In 
this paper we will test this idea on available data and use the results to predict shadowing 
for protons. We do not try to disentangle the two sources of nuclear shadowing, mentioned 
above, gluon recombination and multiple scattering. 

We begin by testing the universality of the relation between the number of gluon partic- 
ipating in the interaction and the ratio of the nucleus and nucleon structure functions. The 



general expression for the nuclear photoabsorption cross section 0, [12|, [18[ reads 

a(p,x)T(b)' 



at A (x,Q 2 ) = 2jd 2 b/l- 



1 



2A 



(2) 



The averaging here is the same as defined in eq. ([I]). T(b) ~ J^dzpA^b, z) is the nuclear 
thickness function, where pA{b, z) is the nuclear density, which dependents upon the impact 
parameter b and the longitudinal coordinate z. The cross section for the interaction of a qq 
fluctuation with a nucleon in eq. (0) can according to standard Regge phenomenology be 
represented in the form, a(p,x) = a(p)x~ Ap ^ 2 \ 

Expanding eq. (0) one can represent nuclear shadowing effects in the form, 



Ra/n{x,Q ) = A _^ Nf __ = 1- 7 7^77^ (T(6)) I-) + ... , (3) 



where (T(6)> = (A - 1) /A 2 J d 2 bT 2 (b). 
Assuming cr(p) ~ C p 2 we obtain 

<a 2 (p)> _ 2AC 



(4) 



(<r(p)) 2m2/n(g 2 /2m2) 
Eq. (Q) looks very similar to the well known Glauber approximation |2(J . It is, however, 
very different. In the Glauber approximation the first rescattering correction is proportional 
to (<y(p)) oc 1/Q 2 - In our case it is proportional to eq. (f|) and is of the order of 1/m 2 . This 
comparison demonstrates that, in terms of multiple scattering theory, nuclear shadowing in 



DIS is dominated by inelastic shadowing while the Glauber eikonal contribution [pO 
vanishes at high Q 2 JIB] . 

An important ingredient of eq. (0) is the assumption that the lifetime of a photon 
fluctuation in the nuclear rest frame is long enough that it can be regarded as propagating 
through the whole nucleus with a frozen intrinsic separation p. This is the same as saying 
that x is sufficiently small, x 1/i71nR a , tha all the parton clouds of nucleons with the 
same impact parameter overlap in the infinite momentum frame 0. However most available 
data are in the transition region, where the lifetime, usually called the coherence time, is 
comparable to the nuclear radius. This can be taken into account by introducing a phase 
shift between qq wave packets produced at different longitudinal coordinates, in the same 
way as was done in refs. P2] , [T1J . This is equivalent to the replacement of the mean nuclear 
thickness function in eq. (H) by an effective one 



(f(b)) = ^H Z b 



dz pA(b, z) e 



tqz 



(T(b)) F*(q) . (5) 



We use a Gaussian form for the nuclear density in order to make a factorization of 
expression @ possible. Generally we use the realistic nuclear density of ref. f23fl . FA{q) = 
exp(—q 2 R A /6) is the nuclear longitudinal formfactor, where R\ is the mean square nuclear 
radius. The decrease of the effective nuclear thickness function at large q can be interpreted 
as a result of the short path length of the hadronic fluctuation in the nucleus if we sit in the 
latter's rest frame, or as an incomplete overlap of the gluon clouds of the nucleons which 
have the same impact parameter in the infinite momentum frame of the nucleus. 

In order to calculate the longitudinal momentum transfer in DIS, q = (Q 2 + M 2 )/2v, 
one needs to know the effective mass of the produced qq wave packet. However a qq state 
with definite separation p does not have a definite mass. We evaluate q ~ 2xm^ assuming 
M 2 ~ Q 2 . Thus X IS cL parameter which controls the value of (T(b)). 

According to eqs. (|3[) and (Q) the number of gluons which interact with the qq fluctuation 
during its lifetime is 

1 (ct 2 (o)) /1\ Ap(Q 2 ) 

n(x, Q 2 , A) = lJ^yj(T(b))F 2 A (q) (-) . (6) 

We expect nuclear shadowing to scale as a function of this new variable n(x, Q 2 , A). Now 
we are in position to test this prediction. 

First of all we fixed Ap(Q 2 ) by a fit to recent data from the ZEUS ||10[ and HI [[II] 



experiments on the proton structure functions using a parameterization x~ Ap at fixed Q 2 . 
The results of the fit are plotted as circles (HI) and squares (ZEUS) versus Q 2 in Fig. 1. 
We also included a hadronic point which we fixed conventionally at Ap(0) = 0.1 which 
corresponds to results from various fits of the energy dependence of hadronic cross sections. 

The Q 2 -dependence of the Ap(Q 2 ) as shown in Fig. 1 was fitted by the ansatz ctp(Q 2 ) = 
a + bexp(—cQ 2 ). The solid curve corresponds to a = 0.36 ± 0.016, b = —0.26 ± 0.016 and 
c= 0.052 ±0.009. 

The expected universal dependence of the nuclear shadowing on n(x, Q 2 , A) is not af- 
fected by the overall normalization of this variable. Provided that Q 2 S> m 2 one can test 



this scaling unambiguously despite the uncertainty in the value of the factor C and the 
quark mass m q displayed in eq. (|). 

The values of the variable n(x, Q 2 ) have been calculated from eq. using the data from 
the NMC [p], and E665 experiments, as well as the results of our fit to Ap(Q 2 ). The 
data on the ratio of the nuclear and nucleonic photoabsorption cross sections Ra/n{x,Q 2 ) 
is plotted against the new variable n(x, Q 2 , A) in Fig. 2. The data demonstrate an excellent 
scaling in n(x, Q 2 , A). The concrete choice of the parameters C and m q (which are strongly 
correlated) in eq. was adjusted to fit eq. (H), the solid curve, to the data. We fixed 
2m 2 g = 0.1 GeV 2 and C /2m 2 q = 22. 

We should comment on this procedure: 
(i) Our considerations are valid for small x, so we limit the x-region to x < 0.07. At this x 
the nuclear structure function shows a small enhancement relative to the proton one, what 
results in Ra/n(x, Q 2 ) > 1 for rafiO. For this reason we renormalized eq. (3) by 3%. 

(ii) Data points 0, |J for Q 2 < 0.5 GeV 2 , were excluded from the analysis because they 
are in the realm of the vector dominance model, rather than DIS. They should correspond 
to the same nuclear shadowing experienced by a p-meson as in the real photon limit. This 
is the reason for the saturation of nuclear shadowing at small x, claimed in P, |2fl. On the 
contrary, nuclear shadowing in DIS at small x, but high Q 2 is not supposed to saturate. 
Ra/n(x, Q 2 ) is expected to decrease logarithmically down to small x below the real photon 
limit due to gluon fusion [R, IIR [TJ 



(iii) The data in Fig. 2 show that Ra/n(x,Q 2 ) depends to a good accuracy linearly 
upon n(x, Q 2 , A) for n < 0.3. On the one hand, this implies that higher order terms in the 
expansion in eq. (H) are small. On the other hand it confirms the validity of the small-p 2 
approximation for a(p) used in eq. ([!]) and eq. @. Data for heavier nuclei and/or smaller 
x will probably demand higher order terms in eq. @. 

The shadowing described by eq. (f|) depends entirely on the gluon thickness function of 
the target. It should, therefore, be possible to apply the same formalism to the shadowing 
in DIS off protons. The nuclear shadowing fixes the parameter C/m 2 of eq. (f|), and thus 
it reduces ambiguities faced by our evaluation of the unitarity corrections to the proton 



structure function. For numerical calculations we used the mean thickness function of 
proton (T)Ql/AnB. The value of the slope parameter B « 6 GeV~ 2 was estimated using 
data on pp elastic scattering, assuming the Pomeron factorization. 

The predicted ratio of the measured proton structure function and the bare Pomeron 
contribution, F$(x, (5 2 )/[Ff(x, Q 2 )]p, is presented in Fig. 3 as function of x versus Q 2 . For 
each value of Q 2 we have put constraints on n < 0.3 where nuclear data exist and we are 
sure that the linear approximation @ is valid. 

After subtraction of the unitarity corrections to the proton structure function we can 
determine the bare Pomeron intercept. We refitted the data [III [T(J and deduce the new 
values of A° P (Q 2 ) . The new interpolating dashed curve corresponds to the parameters 
defined above: a = 0.45 ± 0.017, b = -0.35 ± 0.017 and c = .0525 ± 0.008. The function 
A° P (Q 2 ) which we obtained is quite close at large Q 2 to the intercept of the BFLK Pomeron 

Another assumption, which we quietly made above, should be commented upon. As was 
mentioned above, there are two contributions to the shadowing in DIS, one comes from the 
suppression of the gluon density as a consequence of gluon fusion ggQg, which corresponds 
to the triple Pomeron graph in the framework of standard Regge phenomenology. Another 
contribution to the shadowing comes from the Glauber-like rescattering of the qq fluctuation 
off gluons. This process can also be represented as a parton fusion, but as a fusion of gluons 
into a qq pair, ggQqq. In Regge-model language this process corresponds to the Pomeron- 
Pomeron-Reggeon graph. Both mechanisms show the same x-dependence oc 1/ x 2Ap , except 
for logarithmic corrections. The logarithmic corrections for the triple- Pomeron contribution 
differ for nuclei and the nucleon. Our estimate shows, however, that, in the x and A domain 
investigated here, this difference is small and it decreases with xBO. 

To summarize: we have found a new variable n(x, Q 2 ,A) which all available data on 
nuclear shadowing in DIS scale with at small-x. This variable gives a measure of the 
number of gluons which a hadronic fluctuation of the virtual photon interacts with during 
its lifetime. 

The data on nuclear shadowing allowed us to fix the important parameter C/m 2 which 



determines the screening of the virtual photon, and this reduced the uncertainties in the 
estimations of the unitarity corrections to the proton structure function. A value of the bare 
Pomeron intercept was evaluated from the data on the proton structure function corrected 
for shadowing. This intercept is stiil an effective one, valid in the x region reached by HERA. 
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Figure capture 

Fig. 1 The Pomeron intercept in DIS off a proton, ap(Q 2 ) = 1 + Ap(Q 2 ). The data 
points are results of a fit to the data on x-dependence of the proton structure function 



measured by the HI |TTJ and ZEUS [10] experiments at HERA. The solid curve is the 
interpolation of the results of the fit to the effective intercept. The dashed curve corresponds 
to the bare Pomeron contribution after unitarity corrections are extracted from the data. 

Fig. 2 Data on nuclear shadowing at small x from the NMC Jl], § an d E665 || experi- 
ments versus the scaling variable n(x, Q 2 , A) defined in eq. (|6]). The straight line corresponds 
to F - 2 A /F° = 1 - n(x, Q 2 , A). 

Fig. 3 The unitarity/shadowing correction to the bare Pomeron contribution to the 
proton structure function. Each curve corresponds to the straight line in Fig. 2, plotted as 
a function of x versus Q 2 . 



